RELATIONS FOR BERNOULLI-BARNES NUMBERS AND BARNES ZETA FUNCTIONS 
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Abstract. The Barnes ^-function is 



[x + miai H \-m„a„) 



defined for Re(x) > and Re(z) > n and continued meromorphically to C. Specialized at negative integers —k, 
the Barnes f -function gives 

wfiere B^. (x; a) is a Bernoulli-Barnes polynomial, which can be also defined through a generating function that 
has a slightly more general form than that for Bernoulli polynomials. Specializing 6^.(0; a) gives the Bernoulli- 
Barnes numbers. We exhibit relations among Barnes i^-functions, Bernoulli-Barnes numbers and polynomials, 
which generalize various identities of Agoh, Apostol, Dilcher, and Euler. 



1. Introduction 

We define, as usual, the Bernoulli numbers Bk through the generating function 



e 1 K.. 

A fundamental relation of Bernoulli numbers, known at least since Euler's time, is (for n>\) 



(2) £ . ]BjB„^j = -nBn-i - in-\)Bn. 



Much more recently, multinomial generalizations of dD were discovered by Agoh and Dilcher ||Tl|6l. They 

can be \ 
through 



(n) 

can be viewed as relations between Bernoulli numbers and Bernoulli numbers B\ of order n, defined 



Date: 29 January 2013. 

Mathematics Subject Classification. Primary 11B68; secondary 1IF20, 1IM32. 
Key words and phrases. Bernoulli-Barnes number, Bernoulli polynomial, Barnes zeta function. 

This research project was initiated while M. Beck visited the Universite d'Evry Val d'Essonne; he thanks their Laboratoire 
analyse et probabilites for their hospitality. A. Bayad was partially supported by the FDIR of the Universite d'Evry Val d'Essonne; 
M. Beck was partially supported by the US National Science Foundation (DMS-1 162638). 

1 



2 ABDELMEJID BAYAD AND MATTHIAS BECK 

Dilcher and others proved generalized formulas relating with Bj. The first few are ||6] p. 32]: 

^ = -kBk- i-{k-l)Bu (for any k>l), 

(4) B^l^ =k{k-\)Bk_2 + \K^-2)Bk-i + \{k-l){k-2)Bk (for any A: > 2), 

B^^ = -6k{k-l){k- 2)Bk-3 -\lk{k-l){k- 3)Bk-2 -6k{k-2){k- 3)Bk- 1 (for any k>3). 

Our first goal is to derive relations among Bernoulli-Barnes numbers B/t(a), defined for a fixed vector 
a = (^1,(32, •••,««) £ lR>o through 

Note that, with ^ and (l5]l, the Bernoulli-Barnes numbers and Bernoulli numbers are related as 

miH ^in„=k \ ' "/ 

Of course, one retrieves the Bernoulli numbers of order n with the special case ai = = • • • = a„ = 1, and 
the Bernoulli numbers by further specializing n= I. Our first main result is as follows. 



Theorem 1. For n>3, a = (ai,a2, • • • S lR>o> '^'^'^ "JJ k>n — 2, 

^feV J-2 j(^-n + ;)!|,^/'-"+^'^ '^"\0 otherwise, 

where the inner sum is over all subsets I CI {1,2, . . . ,«} of cardinality j, and a/ := (a,- : / G /). 

Corollary 2. For n>3 and odd k>n — 2, 

" fn + j-4\ k\ (n\ n) ^f-l ifn = k = 3, 
j^2\ i-2 ){k-n + j)\\j) |0 otherwise, 

For example, for « = 3,4 and odd ^ > 1 Theorem [T] gives the relations 

= -K_\ (^>3), 

Bf^ = -2kB^^\-k{k-l)Bf\ {k>4). 

More generally, for any positive integer n>4, Corollary |2] gives the following recurrence formula for the 
numbers : 

(6) pn-4\ („)^_'-^Y" + ^'-^^ ^' ^ ^. foranyodd^>.-2. 

The novelty of these relations as, e.g., compared with dD is that they are between Bernoulli numbers of order 
higher than 1 . We suspect that there are relations analogous to (O for even k but leave the search for them 
as an open problem. 

One of the significances of Bernoulli numbers lies in the fact that they are essentially evaluations of the 
Niemann (^-function l^{z) := Lm>i"^^ (meromorphically continued to C) at negative integers —k: 



RELATIONS FOR BERNOULLI-BARNES NUMBERS AND BARNES ZETA FUNCTIONS 3 

Bernoulli-Barnes numbers appear in a similar fashion in relation with the Barnes ^-function 

C„(z,x;a) := £ — ■ 

mti-'io (x + mai-l hm„a„) 

defined for Re(x) > and Re(z) > n and continued meromorphically to C II7J[8j|9j[14l[l7j[T9l. Specialized 
at negative integers —k, the Barnes ^-function gives 

(-l)"/ti 

(7) i;n{-k,x;a) = -———Bk+n{x;a) 

[k + n)\ 

where Bk{x;a) is a. Bernoulli-Barnes polynomial, defined through 

Thus the Bernoulli-Barnes numbers are the special evaluations Bk{a) = Bk{0;a). It is clear that the Barnes 
zeta function is a multidimensional generalization of various Riemann-Hurwitz zetas functions; e.g., when 
n = I and a = (a), the function l^{s;x,a) is the classical Hurwitz zeta function a^%{s;^). Likewise, the 
Bernoulli-Barnes numbers and Bernoulli-Barnes polynomials extend the (generalized) Bernoulli numbers 
and polynomials to higher dimensions. Further generalizations of Bernoulli numbers and polynomials in- 
clude UTTllTSl. 

Our second main result expresses the Barnes zeta function in terms of Bernoulli-Barnes polynomials, 
Hurwitz zeta functions, and Fourier-Dedekind sums {i}, defined as 

ar[ai,. . . ,aj . . . ,aj;aj) :- 



a 



J x". 



Fourier-Dedekind sums generalize and unify many variants of (generalized) Dedekind sums; see, e.g., [ 161 
or El Chapter 8]. 

Theorem 3. Let ai , . . . ,a„ be pairwise coprime positive integers. Then 

as;x, a) = ~ ^) B„-i-,(x; a) - k;x) 

" f x + r\ 

7=1 r=0 V / 

Theorem[3]has many applications. Specializing s at negative integers gives, with the help of Q: 
Corollary 4. Let a^, . . . ,an be pairwise coprime positive integers. Then 



y\af O-riai, . . . ,aj, . . . ,an',aj)Bm+i 



l (m+l)! ,/n-l\ g.+.+ l(x) 

This is reminiscent of a reciprocity law for generalized Dedekind sums, due to Apostol ||2l; this can be 
illustrated more easily in the case n = 2 and a = {a,b), for which Theorem |3] specializes to: 
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Corollary 5. Let a,b be coprime positive integers. Then 
i;{s;x,{a,b)) = -^^{s-l-x) + (l " ^) 

Again the specialization of s at negative integers gives, for « = 2 and a = {a,b), using ([7]l: 
Corollary 6. Let a, b be coprime positive integers. Then 



1 Im+l fx , 

Bm+2{X, {a,b)) + — — —Bm+2{x)+ ]B,„+i[x). 



m + 2 abm + 2 \ab 

This is a "polynomial generalization" of Apostol's reciprocity law m 

-{a'"-'sm{a,b)+b'^-'s^{b,a))= . ^'"+; + . ]. (aB-bBr+\ 

' (m+l)ab m[m+l)ab 

Here m is a positive integer, a and b are coprime, we use the umbral notation 
{aB-bBr+'=Y. 

and 

(9) S„(^J,^) = £|-^|Bm(^ 

are the Apostol-Dedekind sums. The classical Dedekind sums ll5l[T6| are captured by the special case m = 1. 
Thus in some sense, our study can be viewed as a bridge between Euler-type identities and Dedekind-type 
reciprocity laws. 

Finally, we discuss the special case a = (l,...,l)of Theorem |3] Denote 

C,n{s;x) = C,{S\X,{\,...,\)), 

the Hurwitz zeta function of order n. Since in this case the sums o,-{a\, . . . ,aj, . . . ,an',aj) vanish, we obtain 
the following identity. 

Corollary 7. For any positive integer n, 

Specializing once more s = —m at a negative integer gives: 
Corollary 8. For any positive integers n,m, 



m + k+\ 
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above corollaries the coefficients of the polynomials S^'^j^(x) can be explicitly given by Stirling numbers 



Corollary [8] recovers once more Euler's identity (0 and Dilcher's results in ||6]. We also note that in the 
)ove corollaries the coefficients c 
of the first kind s{n,k) as follows: 

(.0) = "t'("';*)»(„,.+*+i)y. 

(see, e.g., HOl Equation (52.2.21)]). 

2. Proof OF Theorem [T] 

Our proof is based on identities of generating functions. Fix a = {ai,a2, ■ ■ ■ ,an) E M" g define for a 
subset / C {1,2, . . . ,«} 

Theorem [T] will follow from a relation of the even/odd parts of the functions F^^\z), and to this extent we 
define 

even part of F^j'^ (z) if j is odd, 
odd part of F^ j\z) if j is even. 

Note that the even/odd part of F^^^ (z) has the compact description 

z^'\ {e'-^'e,a- ± 1) 
|,^.2n,-e/(^''''^-l) ■ 

Propositio„9. |:(";^;^)(-zr^^oo)(.)^{-^^^ l:;^ 

Proof. The case n = 3 is easily verified, so suppose n>4. What the proposition is claiming in this case is 
that the function 

''(^) - i (" )^ 2 ') (-^"-'^"''<^'> = L (" (-z)"-H(z) 

is even, so that if suffices to prove that F (z) equals 



|/|>2 V l-'l ^ / 



|/|>2 

i.e., that 

\I\>2 

is zero. Written with the denominator HiLi {^"'^' ~ 1)> the function F{z) —F{—z) has the numerator 



[\>2 V 1^1 - 2 y V / ..J 



\i 

and so we can rephrase our goal to proving that 



|/|>2 V I-* I ^ / ^ idl 
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is zero. With Yli^i i^"'^' ~ 1) = L/c/C"!)" we can further rephrase our goal to proving that 

(11) I ("l/l'i;') = I £(_i)«-|/H^I,.i:..«, 

|/|>2 V Kl ^ / yc7 |/|>2 V 1^1 ^ / 7c7 

We will show that the coefficients of e^^^'f'^', for any C {1,2, . . . ,«}, on both sides of (fTTI ) are equal. This 
coefficient is on the left-hand side of ([TTI ) equal to 



'S:<-'K'^')Cr^-r) 



The corresponding coefficient on the right-hand side of (fTTT i is 

n-\K\ 



E<-')"-'''-''{7;-r)^s(-')'-'i";^r)("r') 



Thus (fTTI) is equivalent to 



EMK^)(":!^!-o^|:<-)-f;-r)("r)' 

But both sides equal i^^Z^ , as one can prove, e.g., by putting either side into a generating function for k. □ 
Proof of TheoremUl Recalling that a/ = (a,- : / e /), we see that 

•^'^'^ " n,-e/(^"'^-i) " n,-e/(^-"""-i) 

is the generating function for {—l)'^Bk{aj), and we obtain for odd j 



\i\=i 



and for even j 



(-z)"-^£0W(z) = (-1)" £ X (-l)'s^(a/ 



A:' 

' *+n-./odd 

Theorem [Ufollows by extracting the generating-function coefficients of the identity of Proposition |9] □ 

We should remark that Theorem [T] was in part motivated by llT3l in which Katayama proposed a three- 
term generalization of the reciprocity theorem for Dedekind-Apostol sums IH; Apostol's theorem was a 
byproduct of another paper of Katayama llT2l . Unfortunately, the main theorem of llT3l is wrong; to make the 
central integral of the paper work, one has to use the integrand ^^.^^^ .^^^ ^.^^ ^-^ which, unfortunately, 

does not give rise to Dedekind-Apostol sums. However, using this integrand we discovered Theorem [T] 
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3. Proof of Theorem [3] 

The function 

PA{t) ■=#{{ku...,k„) GZ>o: ^ifliH \-knan = t}, 

which counts all partitions of t with parts in the finite set A := {£?!,...,«„}, is called a restricted partition 
function. For example, basic combinatorics gives 

/'{W}(0=^ „-l 

and a slightly less trivial example was proved by Popoviciu [15] (see also [4. Chapter 1]): for a and b 
coprime, 

where {x} =x — [x] denotes the fractional part of x, a^^ is computed mod b, and b^^ mod a. 

The following theorem was proved in [3 |; however, the authors of that paper did not realize the explicit 
role of Bernoulli-Barnes polynomials. 

Theorem 10. Ifai,...,an are pairwise coprime positive integers, then 

P^'^^^ ^ ( _,-., ^n-i{-t;{a\,...,an)) + Y^ O t 
Proof. We give an outline of the proof. As in |,3J, we compute the residues of 



\b-h] 




\a-h\ 


I a J 


[-1 





The residue at z = gives pA{t), whereas the residue at z = 1 gives ^n-\ {—t', (ai , . . . Finally, if 
A is a nontrivial fl,th root of unity, 

Res(F^(z);z = A) ^ 



where the product runs over all 7 = 1 , . . . , n except j = i. Thus 

^ Res(F,(z);z = A) = -a_f (ai, .. ,a„;a;) . 

The residue theorem completes the proof of Theorem [TO] □ 
Proof of Theorem\3l Set t = m\a\ -\ h We can rewrite the Barnes zeta function as follows: 

C(^;^>a)= L 7 ^- r = E7^- 

mi,..^„>0 {x + miai + ---+ mnanY [x + ty 

By applying Taylor's Theorem to the function 1 1-> f;a) att = —x, 

n-l _ A 

B„_i(-f;a) = £(-l)M )Bn-i-k{x;a)ix + tf, 
and so with Theorem [TO] we obtain 

/yi—\\ « 

= ( _ ni I t ]Bn-i-k{x;ii){x + tf + Y,(y-t{ai,...,aj,...,an\aj). 

\^ /t=o V / ;=o 
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Hence 



„\ (-1)" ^VV ( \n 1 . , ^ ^r^ (^-t{a\,.-.,aj,...,a„;aj) 

Q{s;x,a) = _ 2^(-l) B„_i_^(x;a)C(5-^;x)+ 2^ 2- ■ 

V" '^)- k=Q \ / 7=1 f>0 



j=u>o i^+ty 

Note that a^, {a^ , . . . ,aj , . . . ,an;a j) depends only t mod aj. S etting t = maj + r,0<r<aj — \,we obtain 

\n— 1 1 



+ ^ ^ a_^(ai,...,ay,...,a„;a^-) 



Writing (;c + r + maj)^'^ = aj^ ^ ^ + completes the proof of Theorem |3] □ 

4. The special case a = (a, 1, 1, . . . , 1) 

In the special case A = {a, 1, 1, . . . , 1} (with n I's), most of the terms in Theorem [TO] disappear and we 
obtain 

(12) P{aXK..,i}{t) = ^— (a, 1, 1, . . . , 1)) + a_, (1, 1, . . . , l;a) . 

On the other hand, we can apply (4] Theorem 8.8] to this special case; thus for f = 1 , 2, . . . , a + « — 1, 

a,(l,l,...,l;a) = -^^^^^ B„(?;(a,l,l,...,l)). 

Since a? ( 1 , 1 , . . . , 1 ; a) only depends on t mod a, this range for t is enough to determine af ( 1 , 1 , . . . , 1 ; a) : 

fm^B„(f moda; (a,l,l,...,l)) iff^Omodc?, 

(13) a, ( 1, 1, . . . , l;a) = < , ,{„_i 

[i4^B„(a;(a,l,l,...,l)) if? = Omoda. 

For the case ? = mod a we can also use |4, Theorem 8.4] which gives 

ao(l,l,...,l;a) = l-^^B„(0;(a,l,l,...,l)). 

(An easy way to see that our two formulations of ao (1, 1, • • ■ , l;a) are equivalent is through the difference 
formula 

Bm{x + aQ;{aQ,a\,...an))-B,n{x; (c?o,fl:i, . • .a,,)) =mB,„-\ (x; (ai ,^2, . . 
and then specializing this to m = n, x = 0, qq = a, and = a2 = ■ • ■ = an = 1-) Substituting this back into 
(fT2l) gives, with Xa{t) '■= 1 if a\t and Xa{t) '■= otherwise: 

Proposition 11. = (a, 1, 1, . . . , 1)) -B„(? mod a; (a, 1, 1, . . . , 1))) +Xa(0 ■ 

Using ([T3] ) and Theorem [3] we obtain: 
Proposition 12. Let a = (a, 1, 1, . . . , 1), where a is a positive integer. Then 

Specializing s = —m at negative integers gives, by Proposition [12] with the help of (Q, the following 
formula. 
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Corollary 13. Let a = (a, 1, 1, . . . , 1), where a is a positive integer. Then 

5. Difference, symmetry and recurrence formulas FORB„,(A;;a) 
We conclude by giving various formulas for B„,(x;a), starting with the following difference formula. 

Theorem 14. For a = (ai , . . . , a„) € M" q, we have the difference formula 



(-i)"B,„(-.;.)-B..(.:;.) = ».!i; i 

with Bm{x,aj) = x'" if I = 0. Furthermore, 

n 

B„,(^ + ^ai;a) = (-l)'"B,„(-^;a) . 

i=\ 

Proof. By use of the identity 

/C{1,.-,m} iS/ 

we obtain the formula 



(14) 



1=1 



where / = { 1 , . . . , On the other hand, we have the equality 



-tYe" 



(15) = ^r— • 

i=l i=\ 

Therefore, by and ([T5] ). 



1=1 



Together with ([8]) this completes the proof. □ 

Our next result is a symmetry formula. 

Theorem 15. Let a = (ai , . . . , a,,) G M" q with A\= a\-\ h > 0. Then for any positive integers l,m> 1, 

we have 
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and 

(17) , (-1)' 



= (-l)'"+iB/+,„+i(x;a) + (-l)'+iB„+„,+i(-x;a) . 
Proof. Observe that from dU we obtain 



^B/+i(;c;a) = (/ + l)B/(x;a) 



and by applying the operator ^ to ([T6l ). this implies ([TT] ). 
Now we prove ([T6l ). Consider the generating function 

m>0/>0 \ yt=0 V^/ / 



m! Z! 



££A-^B,,,(x;a)^£(-l) 

/>Oyt>0 (>0/t=0 V' 

/>Oyt>0 (>0/t=() W 



Similarly the generating function is also equal to 



e"- — e 



m>0/>0 \ /t=0 V"^/ / ■Q (^^o/z _ ^a/.v 



Specializing Theorem [T5] at 1 = mwe obtain a recurrence formula for the polynomials 

P„,(x) := {m+l)A-'"(B,„{-x;a) +B^(x;a) 

as follows. 

Corollary 16. For any positive integer m>\, 

2m+l -/m+l\ 



In the case n = \, a\ = \, the polynomials Pm{^) <^re reduced to (m + 1) yB,n{—x) +Bm(x)j . 
Moreover, for the Bernoulli-Barnes numbers we obtain the following recurrence formula 
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Theorem 17. For any positive integer m>\, 



and 



1 /m _|_ 1 \ 

ito|/H;t(2'« + l-« + ^)! 

Note that for « = 1, the above results speciaHze to the well-known difference, symmetry and recurrence 
concerning the ordinary Bernoulli numbers and polynomials. 

References 

1. Takashi Agoh and Karl Ditcher, Higher-order recurrences for Bernoulli numbers, J. Number Theory 129 (2009), no. 8, 1837- 
1847. 

2. Tom M. Apostol, Generalized Dedekind sums and transformation formulae of certain Lambert series, Duke Math. J. 17 (1950), 
147-157. 

3. Matthias Beck, Ricardo Diaz, and Sinai Robins, The Frobenius problem, rational polytopes, and Fourier-Dedekind sums, J. 
Number Theory 96 (2002), no. 1, 1-21, arXiv: math. NT/0204035. 

4. Matthias Beck and Sinai Robins, Computing the continuous discretely: Integer-point enumeration in polyhedra. Undergraduate 
Texts in Mathematics, Springer, New York, 2007, Electronically available at http : / /math . sf su . edu/beck/ ccd . html. 

5. Richard Dedeldnd, Erlauterungen zu den Fragmenten xxviii. Collected Works of Bemhard Riemann, Dover Publ., New York, 
1953, pp. 466-478. 

6. Karl Dilcher, Sums of products of Bernoulli numbers, J. Number Theory 60 (1996), no. 1, Ti-A\. 

1. Emilio Elizalde, Some analytic continuations of the Barnes zeta function in two and higher dimensions, Appl. Math. Comput. 
187 (2007), no. 1, 141-152. 

8. Chelo Ferreira and Jose L. Lopez, Asymptotic expansions of the double zeta function, J. Math. Anal. Appl. 274 (2002), no. 1, 
134-158. 

9. Eduardo Friedman and Simon N. M. Ruijsenaars, Shintani-Bames zeta and gamma functions. Adv. Math. 187 (2004), no. 2, 

362-395. 

10. Eldon J. Hansen, A table of series and products, Prentice-Hall, Englewood Cliffs, N.J., 1975. 

11. Ken Kamano, Sums of products of hypergeometric Bernoulli numbers, J. Number Theory 130 (2010), no. 10, 2259-2271. 

12. Koji Katayama, Barnes' double zeta function, the Dedekind sum and Ramanujan's formula, Tokyo J. Math. 27 (2004), no. 1, 
41-56. 

1 3 . , Barnes ' multiple zeta function and Apostol 's generalized Dedekind sum, Tokyo J. Math. 27 (2004), no. 1 , 57-74. 

14. Kohji Matsumoto, Asymptotic expansions of double zeta-functions of Barnes, ofShintani, and Eisenstein series, Nagoya Math. 
J. 172 (2003), 59-102. 

15. Tiberiu Popoviciu, Asupra unci probleme de patitie a numerelor, Acad. Republicii Populare Romane, Filiala Cluj, Studii si 
cercetari stiintifice 4 (1953), 7-58. 

16. Hans Rademacher and Emil Grosswald, Dedekind Sums, The Mathematical Association of America, Washington, D.C., 1972. 

17. Simon N. M. Ruijsenaars, On Barnes' multiple zeta and gamma functions. Adv. Math. 156 (2000), no. 1, 107-132. 

18. Junya Satoh, Sums of products of two q-Bemoulli numbers, J. Number Theory 74 (1999), no. 2, 173-180. 

19. Mauro Spreafico, On the Barnes double zeta and Gamma functions, J. Number Theory 129 (2009), no. 9, 2035-2063. 

Laboratoire analyse et probabilites, Universite d'Evry Val d'Essonne, 23 Bd. de France,, 91037 Evry, 
France 

E-mail address: abayad@maths . univ-evry . f r 

Department of Mathematics, San Francisco State University, 1600 Holloway Ave, San Francisco, CA 
94132, U.S.A. 

E-mail address: mattbeck@sfsu.edu 



